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Abstract 

The Bk parameter is discussed in the general context of calculating 
beyond the factorization approximation for hadronic matrix elements. A 
variant of the 1/N C method of Bardeen et al. is used. We present calcula- 
tions within a low energy approach and a calculation within the Nambu- 
Jona-Lasinio model. Matching with the QCD behaviour and dependence 
on non-zero current quark masses is studied. 
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1 Introduction 



The problem of calculating matrix elements of weak decay operators is a rather 
old one. In this letter we will restrict ourselves to the matrix element between 
two on-shell kaon states of the following AS = 2 four-quark operator 

O as =2(x) = L; d (x)L» d (x) (1) 

with L sd {x) = s(a;)7^ d(x) and summation over colours is understood. 

This matrix element is usually parametrized in the form of the S^-parameter 
times the vacuum insertion approximation (VIA) as follows 

(K°\O AS=2 (x)\K ) = pMf 2 K m 2 K (2) 

where fx denotes the fx — » liv coupling (fx = 113 MeV in this normalization) 
and rrix is the K° mass. The /i-scale dependence of Bx reflects the fact that the 
four-quark operator Oas=2 has an anomalous dimension and its matrix element 
runs between the scale where the operator has a hadronic realization and the scale 
where it matches the quark-gluon realization of QCD. This anomalous dimension 
is known and using the renomalization group leads to the definition of the scale 
invariant quantity 

B K = B K (fi)a s (fi) a+ (3) 

with 

at one-loop. Here (5^ is the QCD beta function first coefficient. For three active 
quark flavours and N c = 3 we have a + = —2/9. 

The vacuum insertion approximation was historically the first way this par- 
ticular matrix element was evaluated |l[|. Here by definition we have Bk{^) = 1 
at any scale and we can only obtain an order of magnitude estimate. Next this 
matrix element was related to the AI = 3/2 part of K — > tttt by Donoghue et 
al.[]2| using SU(3) symmetry and PCAC. This leads to a value B K ps 0.37. It 
was then found that this relation has rather large corrections due to SU(3) 
breaking. Then three new analytical approaches appeared, the hadronic duality 
approach ||], QCD sum rules using three-point functions^ and the l/N c ex- 
pansion |J. The lattice QCD method also started producing preliminary results 
around this time. A review of the situation several years ago can be found in 
the proceedings of the Ringberg workshop devoted to this subject [0. All these 
approaches have in common that they try to get a numerical value for the Bk 
parameter and study its dependence on the renormalization scale /i. All of these 
methods have been updated and refined. The hadronic duality update can be 
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found in ||, a QCD sum rule calculation is in || and the l/N c expansion method 
has had the vector meson contribution calculated in a Vector Meson Dominance 



(VMD) model [10]. A review of recent lattice results can be found in |11|. A full 
Chiral Perturbation Theory (CHPT) approach to the problem is unfortunately 
not possible. The data on kaon nonleptonic decays do not allow to determine 
all relevant parameters at next-to-leading order ((9(p 4 )) in the nonleptonic chiral 
Lagrangian|jl2|. A calculation of these parameters within a QCD inspired model 
can be found in [[H| where the determination of the B^ factor is done to 0(p 4 ). 
The leading order result for Bk in the 1/N C expansion is well known 

B k Qm)=B k = 1 (5) 

This result is model independent. However, to go further in the 1/N C expansion 
requires some model dependent assumptions. Different low-energy models are 
then used in variants on the 1/N C met hod |J. An example is the calculation done 
within the QCD-effective action model[Q. In fact the tt + — 7r° mass difference 
has been used to test all of the different 1/N C methods, the original one including 
vector mesons ||15|| , the QCD-effective action model jl(| and the extended Nambu- 
Jona-Lasinio (ENJL) model [17]]. The advantage of the latter model is that it has 



been shown to give a good picture of the parameters in the chiral Lagrangian 
including the SU(3) breaking ones, see and references therein. In fact this 



model accounts for a rather large body of low energy hadronic data, see |19[ for a 
review. The reason to resort to a model of this type is that if we also would like 
to know possible SU(3) breaking effects, we cannot use a simple vector meson 
model because we need to know the SU(3) breaking of couplings of the resonances 
to external currents. In practice a large part of these effects comes from scalar 
resonances which are experimentally rather poorly known. We could then use 
a model to estimate these couplings but it seems to make more sense then to 
calculate the relevant quantity directly within the model. It is also important 
to know the higher order couplings beyond order p 4 since weak matrix elements 
involve an integration over internal momenta. 

In this letter we shall give the results for the next-to-leading corrections to (|5|) 
for the simple NJL model without a vector like coupling. This does not include 
vector mesons so it will not be the final answer but one can already do a study 
of the SU(3) breaking effects since these tend to come mainly from the scalar 
sector. Results on the full model and more details regarding the present work 



will be presented elsewhere [20]. The full model requires inclusion of anomalous 
diagrams and involves about an order of magnitude more terms than the present 
work. 

We will first discuss the method and then use CHPT to derive first results. 
This is essentially the method of M in a slightly different notation. This will also 
allow us to test at low cut-offs the calculation done afterwards. Next we describe 
the NJL model and give the lowest order result in 1/N C . Then we describe the 
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main part of the work. The calculation of the non factorizable part of the AS = 2 
two-point function in the NJL model. Next we discuss how the Fierz-contribution 
is reproduced analytically in this approach and how it affects the final result. We 
will give numerical results for the case of massless quarks, of non-zero degenerate 
quark masses and of different quark masses. In the end we present our main 
conclusions. 

2 The method and definitions 

We calculate here not directly the P^-factor but the AS = 2 two-point function 

G F n AS=2 ( 9 2 ) = iyrf 4 xe^e irAS = 2 (0|T(P ds (0)P ds (x)) |0) 

= ^ 2 |ci 4 a;e^(0|T(P ds (0)P ds (x)r A5=2 ) |0) (6) 

in the presence of strong interactions. We shall use the NJL model for scales 
below or around the spontaneous symmetry breaking scale. Here Gp is the Fermi 
coupling constant, we use P ds {x) = d(x)i r y 5 s(x), with summation over colour 
understood and 

Tas=2 = -G f J d A yO AS=2 (y). (7) 

The reason to calculate this two-point function rather than directly the matrix 
element is that we can now perform the calculation fully in the Euclidean region 
so we do not have the problem of imaginary scalar products. This also allows 
us in principle to obtain an estimate of off-shell effects in the matrix elements. 
This will be important in later work to assess the uncertainty when trying to 
extrapolate from K — > n decays to K — > 2-k. This quantity is also very similar 
to what is used in the lattice and QCD sum rule calculations of B K . 
The AS = 2 operator in (0) can be rewritten as 

T A s=2 = -G F J #yL?{y)lUv) 

= -G F J -^yj d\, J d±x 2 e-^-^Lf( Xl )L» d {x 2 ). (8) 

This allows us to consider this operator as being produced at the Mw scale by the 
exchange of a heavy X AS = 2 boson. We will work in the Euclidean domain 
where all momenta squared are negative. The integral in the modulus of the 
momentum r in (§) is then split into two parts, 

rM w r-p, rMw 

/ d\r\ = / d\r\ + / d\r\ . (9) 
Jo Jo J n 

In principle one should then evaluate both parts separately as was done for the 
7r + — 7T° mass difference in the above quoted references. Here we will do the upper 



3 



part of the integral using the renormalization group. This results in the integral 
being of the same form but multiplied with the Wilson coefficient C(/z), 



AS=2 



(27T) 



(10) 



A strict analysis in l/N c would correspond to set C(/x) = 1 + d([i)as(n) and 
to evaluate the second term using factorization in leading 1/N C . d(fi) is what 
the one-gluon exchange diagrams would give with a lower cut-off /z. We will 
however use the full one-loop Wilson coefficient C(/i) = (as{n) / ces(M w )) a+ . In 
the remainder formulas this factor will be suppressed for simplicity. The lower 
part of the integral is then evaluated using the NJL model. 

To illustrate the procedure let us calculate the relevant quantities in chiral 
perturbation theory. This will only make sense for a scale [i such that CHPT 
can be trusted. The lowest order (LO) contribution comes from the diagram 
in fig. [l|a. The vertices at this order can be obtained from the standard chiral 
Lagrangian of order p 2 . The result is 



n 



AS=2 



'CHPT-LO 



-2B 2 F 4 



(q 2 



m 



K ) 



The parameter B is related to the quark condensate in the chiral limit by (qq) = 
—B Fq and F is the pion decay constant in the chiral limit, F »s 86 MeV. 
At this order, the only nonvanishing contributions are leading in 1/N C , therefore 
when the result in (O) is properly reduced (see pif ) and compared with eq. (0) 
one gets the leading order l/N c result in eq. (|5]). 

At next-to-leading order (NLO) in the chiral expansion more diagrams have to 
be taken into account. The ones that are non-factorizable (NF), i.e., the ones that 
are not reducible when the X boson propagator gets cut, are depicted in fig. |l]b-f 
(plus the last one symmetrized). Here the vector- vector part of the interaction 
contributes in diagrams b and d while the axial-axial part contributes in c,d,e 
and f. Diagram d was not considered in the usual 1/N C approach. We will discuss 
its significance later. We have to identify the cut-off /i uniquely from diagram to 
diagram. We would also like to get the model independent leading l/N c result 
independent of the value of the scale //. We do this by routing the external 
momentum (q) and the loop momentum (r) through the X boson propagator 
(Px — q + t ) and then cutting on the loop momentum in the Euclidean by 

— Now, in order to compare with the QCD realization in terms of quarks 
and gluons identifying fi as the renormalization scale, this then requires q 2 / fi 2 
and Xq CD / fi 2 to be small. In the chiral limit (massless quarks) all integrals can 
be performed analytically and the result is 



n 



~-2(q 2 )cHPT- 



■NF 



BlF 2 
16tt 2 



4/i 2 



+ 1 



(12) 
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• ■ • • ■ • 

(a) (b) 




Figure 1: Chiral Perturbation Theory contributions to IIas=2 (<? 2 )- (a) Lowest 
order, (b-f) Higher order non-factorizable. A circled cross is an insertion of the 
external pseudoscalar current, P sd (x), a dot is a strong interaction vertex and 
the square represents the AS = 2 operator, T& s =2- The full lines are meson lines. 
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The quartic dependence on \i cancels between the different diagrams as required 
by chiral symmetry. Therefore we get for the Bx factor 



CHPT 



1 - 



16tt 2 F 2 



(13) 



The correction is negative. It disagrees somewhat with the result obtained in || 
because there no attempt at identifying the cut-off across different diagrams was 
made. Since we work at leading level in 1/N C in the NLO CHPT corrections we 
have included the relevant singlet (rji) component as well using nonet symmetry. 
The correction in (|T3|) has precisely the right behaviour to cancel partly C{/a) 
which increases with increasing \i. We will defer a numerical discussion till later. 



3 The NJL model 

The Lagrangian of the NJL model is given by 

Cn jl = q [«7 M {dp - ton - ia^ 5 ) - M — s + ip-f 5 ] q 

+ 2g s Y J {f R q i L ){iJ R ) (14) 

with q — {uds),i,j are flavour indices and M. = diag(m u , m d , m s ). The coupling 
gs is given by g$ = 47r 2 Gs/ (N c A^) in terms of the notation used in [|17|, [18| pl|j . 
v^a^s and p are external vector, axial-vector, scalar and pseudoscalar fields. 
These are used to probe the theory with. This Lagrangian can be argued to follow 
from QCD in the following way: all indications are that in the pure glue sector 
there is a mass gap. The lowest glueball in lattice QCD numerical simulations has 
a mass of about 1.5 GeV. This means that correlations below this scale should 
vanish. We then treat the interactions of quarks below this scale as pointlike. It 
can also be seen as the first term in an expansion in local terms after integrating 
out fully the gluons. The Lagrangian in ( P~3| ) has the correct chiral symmetry 
properties. Thus results calculated within this formalism will have the correct 
chiral properties. 

This model, for Gs > 1, spontaneously develops a quark vacuum expectation 
value so that the chiral symmetry is spontaneously broken. We will now work 
within the 1/N C expansion in this model. How to calculate two-point functions 
and three-point functions in the presence of non-zero current masses can be found 
in ||21|| . We have for this work evaluated all necessary two-, three- and four-point 
functions in terms of different masses. The prescription used for calculating 
the integrals is the same as in ET[. We first use identities of the type 2q ■ r = 



((q + r) 2 -M 2 )-q 2 -(r 2 -M 2 ) and r 2 = (r 2 -M 2 ) + M 2 to reduce all integrals to 
scalar loop integrals. In these, we then combine the propagators using Feynman 
parameters. Then we rotate to Euclidean space and regularize the integrals using 
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Figure 2: The leading 1/N C contribution to Has=2{q 2 ) in the NJL model. Symbols 
as in fig. m except that a dot is now a NJL- vertex and the full lines are constituent 
quark-lines. The flavour is mentioned next to the lines. 



a consistent proper time cut-off which introduces a cut-off scale A x in all the loops 
with constituent quarks. This procedure reproduces the two-point functions of 
[|17jj and in a low-energy expansion all results of |18[. The flavour anomaly is also 



correctly treated in this way, see p2 



As input parameters we use those obtained in the best fit to F and the Lj of 



CHPT in |T8[ for the full ENJL model, i.e. including the vector-axial- vector four- 
quark term which is modulated by the coupling Gy. At present we are mainly 
interested in seeing the effect of SU(3) breaking and we shall put Gy to zero, i.e. 
we work within the model in eq. The use of the same input values here as 

the case with Gy 7^ will also allow to get a better comparison with that more 
realistic case which is in preparation. These input values are Gs = 1.216 and 
A x = 1.16 GeV. This gives a constituent quark mass of about 265 MeV in the 
chiral limit. We will use = 3.2 MeV and m s = 83 MeV. These are the values 
that reproduce the experimental pion and kaon mass using the full ENJL model. 

In the NJL model of eq. (fUj), these input parameters correspond to the 
quark condensate (gg)chirai ~ -(281 MeV) 3 , F « 114 MeV, B « 1.7 GeV and 
m 2 K ~ 0.133 GeV 2 . We will also study the case where both quark masses are 
equal because this is where the lattice calculations have been done so far, then 
we shall use m s = = 43 MeV. This leads to the same kaon mass as given 
above. 



4 The NJL calculation 



We now proceed to evaluate 11^5=2 i n the NJL model. The leading contribution 
comes from the class of diagrams in fig. 0. In terms of the mixed pseudoscalar- 



axial-vector two-point function Hp^(q) s d (see |JT7 ] for definition), the result is 

,(g 2 ) = -In^(-g) sd n PM (g) sd . (15) 



n 



AS=2{ 



This mixed two-point function has been calculated in the full ENJL model in the 
chiral limit |T7] and in |?T| for non-zero quark masses to all orders in the chiral 
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Figure 3: The class of 3-point diagrams. Symbols as in fig. |^, the hatched areas 
are a summation over sets of one-loop diagrams as shown in fig. ||. 



expansion. The lowest order term in the expansion of the result in (|15|) coincides 
with that in ( |TTD when the parameters are those determined from the NJL model. 

At the non-factorizable level life becomes a lot more complicated. There are 
two major classes of diagrams in the leading contribution to the non-factorizable 
case. They are depicted in figs. |3| and f|, we will refer to these as 3-point and 
4-point diagrams respectively. In addition the diagram of fig. |3] also exists with 
everywhere the quarks s and d interchanged and the directions of all fermion 
lines reversed. The latter diagram can then related to the first version by CPS- 
symmetry and can be calculated using the same code written for the diagram of 
fig. |3] but with the quark masses m s and interchanged. We have checked that 
diagram in fig. [| gives the same result after that interchange. The tails here can 
be 0,1,2,3,... loops and they can be resummed using the procedure given in [|17|. 
The expressions needed for the tails are given in [^1|] . They represent in a sense the 
fermionic equivalent of the propagators of the mesons. Notice that the two-point 
functions needed also contain parts that do not have a pole. Here it can be seen 
where the main complexity of the calculation comes from. The X boson has both 
vector and axial-vector couplings, the internal coupling of the tail to the 3 or 4 
point quark loop can thus have 4 different Dirac structures. The three- and four- 
point functions themselves have a rather complicated momentum dependence 
and need an explicit numerical integration over the Feynman parameters to be 
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Figure 4: The class of 4-point diagrams, symbols as in fig. [| 



evaluated. This makes it rather important to have several consistency checks 
on the final numerical result. We will describe the ones we have used below in 
addition to the interchange of masses on the 4-point diagram. We refer to the 

3- point diagrams by the form of the insertion at places A,B,C and D and the 

4- point ones by their insertions at A and B. This in turn fixes the insertion in 
the points E and F in both diagrams. As an example let us quote the expression 
for the PSSP 3-point diagram corresponding to the flavours in fig. 

<PSSP> = l -[l + g s U P (q) sd } 2 gl J^-^UU-q-rUUp^q + rU 



(27T) 

The others are considerably more involved. The two- and three-point functions 



x Tlp d p S {q + r, -r) [1 + g s ^-s{-r)dd] H-pps^ - r, r) . (16) 



here are defined with the notation used in |21| . The non-barred ones correspond 
to the full functions and the barred ones to the one loop expressions Hi~7|, |2l| . We 
have made sure that the expressions for the underlying 3 and 4-point one-quark- 
loop expressions obey all the Ward identities that follow from the symmetry in 
the one-loop case for the most general case, all masses different. The current 
identities are satisfied with the constituent masses. At the present level there is 
no contribution from diagrams that produce a Levi-Civita symbol. These only 
contribute at the level considered here as soon as Gy ^ 0. A large fraction of 
the work was involved in making sure the 3- and 4-point one-loop functions were 
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correct. 



5 What happened to the Fierzed term ? 

In the standard VIA we have Bk([i) = (3/4) (1 + 1/N C ). From the diagram in 
fig. H we get B K (ii) = 3/4 and the difference comes from the colour suppressed 
combination of attaching the operator Tas=2 to the quark lines. In terms of the 
X boson propagator this type of contributions corresponds to taking the diagram 
of fig. |] and replacing the lines connecting A-E and B-F by a pointlike insertion 
and sending the cut-off /i to infinity. We can then rearrange the integration over 
r and the internal quark loop momenta such that the new integration variables 
have no mixed dependence in any of the quark propagators. We can now split 
the double integral into two independent parts. The expression for this diagram 
then reduces to the one in fig. ^ with an extra factor of 1/N C in front obtaining 
the usual VIA result. One can wonder how good is the VIA (i.e., the Fierzing) 
when one introduces the cut-off /i in the r momentum. The answer is that this 
l/N c is only reached for rather high values of the cut-off [i. 

Let us also see what happens in the context of chiral perturbation theory. 
In fact the same contribution appears in the CHPT calculation. It is the one 
depicted in fig. [I|d. To the order we have calculated here (i.e. next-to-leading 
order) it is exactly zero, the vector-vector and axial-vector-axial-vector cancel 
exactly. It is only at higher orders in the CHPT expansion, when the vertex 
vector(vector-axial)-V boson is the one coming from the L\q term in the 0(p 4 ) 
chiral Lagrangian, when this contribution starts showing itself. This in fact allows 
us a rare check on a single diagram. The chiral calculation for the vector(axial- 
vector)-vector(axial-vector) piece gives (at NLO) a contribution of 

n a 5 rf fi-5.4,(9 2 ) = (-)|^4^. (17) 

In the NJL calculation, this agrees for small values of fi (up to about 0.35 GeV) 
well. For larger values of \x the contribution of this diagram is (3/4) e, with e ~ 
0.006, 0.03, 0.08, 0.1, 0.2 and 0.3 for fi = 0.3, 0.5, 0.7, 0.9, 1.1 and 1.5 GeV. 
As can be seen the Fierz term is slightly suppressed for the values of fi that are 
relevant in the next section. 

6 Numerical results 

We have studied three cases, namely, the chiral case mentioned above; the case 
with SU(3) symmetry breaking m s ^ with the input parameters introduced 
in section 3 and the case with non-zero quark masses but with m s = also 
introduced in section 3. The procedure we have followed to analyze the numerical 
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H (GeV) 




°K 




B a M 


Dm 
B K 




£>cq 
D K 


0.3 


0.68 


0.50 


0.74 


0.50 


0.55 


0.74 


0.55 


0.5 


0.59 


0.59 


0.71 


-0.44 


0.71 


0.72 


0.72 


0.7 


0.53 


0.58 


0.69 


-2 


0.75 


0.68 


0.75 


0.9 


0.48 


0.55 


0.66 


-3 


0.76 


0.65 


0.75 


1.1 


0.45 


0.54 


0.64 


-4 


0.76 


0.64 


0.76 



Table 1: Results for Bk and Bk- 

results is the following. We fit the ratio between the correction and the leading 
l/N c result for a fixed scale /i to a/q 2 + b + cq 2 which always gives a very good 
fit (a, b and c are /i dependent). Once we have this fit we can extrapolate our 
Bk form factor (remember that we have calculated it for Euclidean q 2 ) to the 
physical Bk, i.e. to q 2 = m 2 K . 

Let us first treat the chiral or massless quarks case. Here a nontrivial check 
on the results is that the diagrams have a behaviour which sums to 1/q 2 . The 
individual contributions do not have this behaviour. In fact for our lowest q 2 the 
cancellations are rather significant. This cancelation happens point by point in 
the loop momentum (r) integral since there is an underlying chiral relation for the 
4-point amplitude. Therefore we have used a deterministic numerical integration 
scheme to have all diagrams evaluated at the same internal momenta. This 
increases the accuracy in the sum beyond the accuracy in the separate diagrams 
since only the noncanceling contributions survive at all stages of the calculation. 
We also see here a small q 2 dependence in the Bk form factor, this is for Bk as 
defined above but with off-shell kaons. 

In this case a is compatible with zero as required by chiral symmetry, b is the 
relevant contribution to Bk since rn 2 K \ x = 0. The first three columns in table [TJ 
are ii, B k {li) and B K for this case. For calculating B K we use a + = —2/9 and 
Ag_ = 250 MeV which corresponds to a^(1.5 GeV) = 0.39. For very small ll 
we also see agreement with the CHPT calculation given above with the values for 
Bq and Fq expected for our NJL input values. The matching with QCD leading 
logarithmic evolution, i.e. the stability of Bk, is not very good in this case. 

In the second case (m s ^ rrid), because the chiral symmetry is broken, there 
is a possibility for contributions to B K (q 2 ) that are not proportional to q 2 , i.e. 
a / 0. In fact a CHPT calculation like the one above predicts precisely the 
presence of this type of terms. For small values of q 2 the part due to a dominates 
even though it is only a small correction when extrapolating to the physical B 7 ^ 
at q 2 = m 2 K - This can be found in column 4. The fifth column is the form factor 
B K for q 2 = —0.001 GeV 2 where the correction due to the a term is sizeable. 
Notice the difference between these two columns. This same feature should be 
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visible in the lattice calculations as soon as they are done with different quark 
masses. The invariant 13% for this case is in column 6. 

In the the last case, i.e. irtd = m s , which is similar to the present lattice QCD 
calculations, the fit gives a = to a good precision and the value of B^ extracted 
is rather independent of q 2 . The invariant B e ^ in this case is in column 8. 

A reasonably stable value for B\ occurs for scales \i ~ (0.7 ~ 1.1) GeV. But 
the difference with the massless case is numerically significant. We obtain 

&k(™>k ~ °- 13 Gey2 ) ~ L35 Bl(m 2 K = 0) (18) 
for scales fi « (0.7 ~ 1.1) GeV. 



7 Conclusions 

We have studied a nonleptonic weak matrix element in a way where the short and 
long distance contributions are separated in a physically visible fashion through 
the X-boson propagator. This allows us to identify the role of \i across diagrams 
and to connect with the running in the short distance domain. We have done 
this in the NJL-model to all orders in momenta there. This gave us a result that 
has manifestly the correct chiral properties. 



In view of the results of |Tq, 19, 21] we expect to get a good prediction for 



the effects of non-zero and different current quark masses. We see those and 
find a significant change due to both. For the extrapolation to the kaon pole the 
difference between the masses has a much smaller effect than the fact that they 
were non-zero. In order to compute Bx in the general case a careful extrapolation 
to the poles was needed. The final correction to the B^ parameter compared to 
its leading value of 3/4 turns out to be rather small. The stability of Bk when 
quark masses are present, is around scales fi ~ (0.7 ~ 1.1) GeV. For these scales 
one may expect that the contribution of spin 1 mesons (remember that only 
the spin mesons are the ones included for Gy = 0) is important. This may 
be implemented like in the case of the tt + — tt° mass difference [15, [16], |17]] by 



going to the full ENJL model with vector-axial- vector mesons included. This also 
will allow us to compare with the results of the next-to- leading 1/N C corrections 
obtained in similar approaches |TIJ. We are studying now this case PD| . We are 



also studying whether we can use the present program to obtain more information 
about higher order terms in the chiral Lagrangian. 
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